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1.  Introduction 

The  interest  to  study  of  boundary  layers  with  blowing  of  reacting  chemical  substances  is 
due  to  extensive  applications  of  this  process.  The  mechanism  of  how  the  porous  blowing  effects 
heat  and  mass  transfer  and  friction  was  studied  deeply  for  nonreactive  flows  [1,  2],  but 
combustion  makes  the  whole  pattern  more  complex.  The  heat  release  in  the  reaction  zone 
creates  a  strong  effect,  which  changes  significantly  the  thermophysical  properties  and. 
consequently,  turbulence  characteristics.  This  was  proved  by  existing  experimental  data  on 
turbulent  structure  and  heat  and  mass  transfer  in  near-wall  turbulent  pulsations  in  the  near-wall 
diffusive  flame  [3-7], 

A  certain  progress  has  been  achieved  in  description  of  this  problem,  and  simulation  of 
turbulent  combustion  in  a  boundary  layer  [8,  9]  gives  analysis  of  different  models  of  chemical 
reactions  [10]. 

However,  this  problem  creates  many  questions  that  still  require  a  more  detail  study.  One 
of  the  issue  is  how  the  intensity  of  fuel  blowing  through  a  porous  surface  affects  the  structure  of 
the  boundary  layer,  position  of  flame  front,  heat  and  mass  transfer,  and  friction. 

The  current  paper  presents  results  of  simulation  for  a  turbulent  flow  with  hydrogen 
blowing  into  the  boundary  layer.  Since  the  process  has  many  parameters,  the  main  emphasis 
was  put  on  analysis  of  dynamic  parameters  -  velocity  profiles,  integral  scales  and  surface 
friction. 


2.  Problem  statement  and  main  equations 


A  plain  porous  plate  for  supply  of  pure  hydrogen  is  exposed  to  a  non-gradient  airflow 
with  velocity  u0  =  20  m/s.  The  hydrogen  fed  through  the  wall  reacts  with  air  oxygen  inside  the 
boundary  layer,  at  the  combustion  front,  and  reaction  products  diffuse  toward  the  external  limit 
of  the  boundary  layer  and  toward  the  wall.  To  exclude  the  temperature  factor,  all  calculations 
were  performed  at  a  fixed  length  and  the  same  temperature  for  the  wall  and  stream  core  T0  =  Tw 
=  300  K.  The  concentration  of  the  blown  hydrogen  was  taken  as  a  constant  over  the  length 
1 )  =  const .  Correspondingly,  the  mass  flux  over  the  plate  length  decreased  with  a  growth 

of  the  boundary  layer  thickness. 

The  intensity  of  hydrogen  feed  was  varied  through  a  change  of  hydrogen  concentration  on 
the  wall.  Our  simulation  was  performed  for  a  most  wide  range  of  concentrations 
[ch2  )  =  0,02  -  0,8.  At  a  lower  concentration,  there  was  no  flame  at  all,  and  we  failed  to 


obtain  a  stable  solution  at  higher  concentrations  of  hydrogen.  This  range  of  concentration 
corresponds  to  a  wide  range  of  the  blow  parameter  bl=2  jw  /  PqUqc  f  =  0, 1  -t-4. 


This  work  uses  k-£  model  of  turbulence  (LRN  -  modification  by  Chien  [12]). 
equation  for  each  variable  can  be  presented  in  the  generalized  form: 
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Variables  <t>,  diffusion  coefficients  T®,  and  source  terms  T*  are  presented  in  the  table. 
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Here  II  -  longitudinal  velocity,  H  -  enthalpy  of  the  mixture,  c,  -  mass  fraction  of  z'-th 
component,  k  -  kinetic  energy  of  turbulence  and  8  -  dissipation  rate  of  turbulence.  Turbulence 
model  constants  and  functions  were  taken  from  [12]. 

The  parameters  in  these  equations  like  viscosity,  Prandtl  number,  Schmidt  number  were 
calculated  by  semiempirical  formulas  for  a  multicomponent  gas  mixture  [11],  and  the  density 
was  calculated  from  the  ideal  gas  law. 

As  preliminary  calculations  demonstrated,  combustion  under  these  conditions  is 
determined  mainly  by  the  rate  of  reactants  mixing,  since  the  reaction  rate  is  rather  high. 
Therefore  we  imply  in  further  calculations  that  only  the  reaction  fuel+oxidizer— ^products  takes 
place,  and  its  rate  is  infinitely  high.  This  corresponds  to  the  so-called  diffusion  model  of 
combustion. 

Under  these  assumptions  the  Lewis  number  is  different  from  one  (Le  1 )  and.  according 
to  data  from  [10],  this  approximation  gives  a  good  agreement  with  experiment  and  the 
combustion  under  these  conditions  can  be  considered  as  a  chemical  equilibrium  process. 

The  values  of  turbulent  numbers  of  Prandtl  and  Schmidt  were  taken  constant  over  the 
boundary  layer  thickness  and  equal  to  0,9. 

The  boundary  conditions  for  this  problem  were  the  following: 

—  on  the  initial  cross-section  Cv  =  0) 

k  =  koj  £  =  £}•  u  =  Uo,'  V  =  Vo,*  T  —  Tq\ 

—  on  the  external  side  of  the  boundary  layer  (y  =  S) 

u  =  uo;  T  =  T0;  C()n  =0.23;  CN^  =0,77  and  PP/chy  lor  other  variables; 


—  on  the  wall  (y  =  0) 


0;  k  =  0;  £=();  C„2  =(CHi)w, 


T  —  Tw;  (jw) 77  — 


pD 
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H  , 


dy 
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The  advantage  of  the  used  turbulence  model  is  that  it  works  at  low  Reynolds  numbers. 
Therefore  it  was  possible  to  calculate  the  flow  even  for  laminar  range.  This  work  uses 
coordinate  transformation  and  numerical  method  for  solving  boundary  layer  equations  that  have 
been  proposed  by  Denny  et.  al.  [13]. 

Before  the  simulation,  we  performed  some  test  runs,  for  checking  of  operability  of 
numerical  algorithm  and  turbulence  model.  The  friction  and  heat  transfer  on  a  impenetrable 
plate  were  calculated  for  blow  of  the  same  (or  foreign)  gas.  Calculations  demonstrated  a  good 
correspondence  with  the  laws  for  friction  and  heat  and  mass  transfer  known  from  literature  [1,2]. 


202 


3.  Simulation  results  and  discussion 

The  friction  coefficient  as  a  function  of  Reynolds  number  Rex  with  hydrogen  blow  and 
combustion  in  the  boundary  layer  is  plotted  in  Fig.  1.  Here  the  curves  for  different  hydrogen 
concentrations  on  the  wall  are  presented,  and,  correspondingly,  data  for  different  intensity  of  blow. 

Let  us  consider  the  behavior  of  the  friction.  With  an  increase  in  the  wall  concentration  and 
blow  intensity,  the  friction  decreases.  However,  unlike  for  nonreactive  flows,  combustion 
creates  a  lag  in  laminar-turbulent  transition:  for  a  low  intensity  of  blow  it  occurs  at  Reynolds 
numbers  Rex  =  (3-f5 )  1 06,  at  this  is  by  one  order  higher  than  for  the  no-friction  situation.  For  a 
higher  intensity  of  blow,  the  transition  shifts  to  the  zone  of  smaller  Reynolds  numbers. 

Existence  of  intensive  heat  release  in  a  boundary  layer  influences  the  distribution  of 
friction  over  the  plate  length.  As  analysis  revealed,  for  laminar  zone  of  flow  the  exponent  in  the 
friction  law  is  different  from  0.5;  for  turbulent  regime  the  blow  creates  variation  in  the  exponent 
in  the  limits  n  =  OvO,  1 6. 

The  profiles  of  velocity,  total  enthalpy,  and  generalized  concentration  are  almost  coincide 
between  them,  and  this  tell  us  about  similarity  for  processes  of  heat  and  mass  transfer  and 
friction  even  under  such  complicated  conditions.  The  prove  that  Reynolds  analogy  is  valid  for 
reacting  flows  (  with  the  mentioned  variables)  was  provided  also  experimentally  [3]. 

Evolution  of  velocity  profiles  over  the  length  of  porous  plate  under  combustion  condition 
is  demonstrated  in  Fig.  2.  Here  the  velocities  are  represented  in  traditional  semi-logarithmic 

coordinates.  Here  u+  =  u I uT\  uT  =  Jtw  /  pw  ;  y+  =  yuT  / v  . 

For  Reynolds  number  Rex  <  106  ,  the  velocity  profile  has  no  the  traditional  intrinsic  zone 
with  a  logarithmic  law.  This  manifests  that  the  flow  mode  is  close  to  laminar;  this  is  caused  by  a 
strong  heat  release  in  the  boundary  layer.  With  a  higher  Reynolds,  the  profile  deviates  from  the 
laminar  distribution,  so  the  flow  becomes  turbulent.  These  tendencies  were  noted  in  Fig.  1  also. 

The  effect  of  blow  intensity  on  velocity  profile  in  a  boundary  layer  is  shown  in  Fig.  3  in 
semi-  logarithmic  coordinates. 

In  this  graph,  the  simulation  results  with/without  combustion  are  compared  directly  (all 
other  external  parameters  are  identical).  At  low  wall  concentrations  (Fig.  3a  -  small  blowing) 
the  velocity  profile  is  close  to  the  standard  distribution  in  a  laminar  sublayer  and  at  the  intrinsic 
logarithmic  zone.  With  combustion,  the  profile  undergoes  a  considerable  change.  As  it  was  for 
low  Reynolds  number,  we  have  no  any  logarithmic  zone:  this  testifies  that  the  influence  of  heat 
release  on  the  flow  structure  is  overwhelming. 


Fig.  1.  Friction  coefficient  in  a  boundary  layer 
with  hydrogen  injection  and  combustion 
(hollow  symbols  -  combustion,  solid  - 
injection  without  combustion). 


y 

Fig.  2.  Velocity  profiles  (combustion  case,  CHi=0.5) 
dependence  on  Reynolds  numbetj.  Solitjl  lines  a^e  stan¬ 
dard  turbulent  ^elocity  profile  u  =  y  ,  if  y  <11, 
u  =  5.75 lg  y  +  5.5  elsewhere. 
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Fig.  3.  Velocity  profiles  in  a  boundary  layer  with  injection  and  combustion  (hollow  symbols  -  combustion, 
solid  -  injection  without  combustion).  Rex  =107 

The  existence  of  flame  front  makes  the  profiles  closer  even  at  high  rates  of  blowing.  We 
can  see  in  Fig.  3b  that  this  difference  became  minimal,  but  the  velocity  profile  (no  combustion) 
display  a  tendency  of  deviation  from  the  classical  logarithmic  distribution. 

There  is  a  strong  effect  of  heat  release  on  flow  structure  at  the  flame  front.  Depending  on 
the  blow  intensity,  the  front  coordinate  may  shift  within  the  boundary  layer  thickness.  This  is 
illustrated  by  data  presented  in  Fig.  4.  With  a  higher  intensity  of  blow,  the  front  coordinate 
moves  out  from  the  wall,  and  this  shift  may  be  considerable.  Obviously,  the  fact  where  the 
flame  front  is  -  in  the  near-wall  or  external  zone,  must  influence  the  transfer  characteristics. 

All  of  these  peculiarities  of  flow  structure  in  a  boundary  layer  with  combustion  can  be  applied 
to  the  behavior  of  the  friction  ratio.  The  dependency  of  this  coefficient  on  blow  parameter  is 
shown  in  Fig.  5.  Flere  C/Cfo  is  the  friction  function,  Cp  is  the  friction  for  a  impenettable 
isothermal  plate,  and  bl  =2  jw  /  PqUqC f  is  the  blow  parameter.  From  data  in  Fig.  5  we  can 

estimate  the  contribution  of  blow  and  combustion  into  reduction  of  friction.  The  influence  of 
every  of  these  factors  reduces  the  friction,  and  this  effect  is  considerable.  With  the  joint  action 
of  both  factors,  the  reduction  in  friction  may  reach  a  high  value.  Flowever,  for  intensive  blow, 
the  results  of  simulation  for  flows  with  and  without  reaction  become  closer  and  the  effect  from 
combustion  is  weaker  because  the  flame  front  is  pushed  out  from  the  wall  and  from  the  outer 
side  of  the  boundary  layer. 
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Fig.  4.  Flame  front  position  dependence  on 
blowing  . 


Fig.  5.  C/Cp  dependence  on  blowing  parameter 
(hollow  symbols  -  combustion,  solid  - 
injection  without  combustion).  Rex=10 
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c) 


Fig.  6.  Comparison  between  experimental  [5]  (dots) 
and  simulation  data  (lines).  x=120mm. 
a)  temperature,  b)  mean  velocity,  c)  turbulent 
pulsations. 


While  comparing  our  simulation  results  with  experimental  data  on  hydrogen  combustion 
[5]  (which  are  qualitatively  similar),  we  faced  the  problem  described  in  [14].  The  matter  is  that 
those  experiments  were  conducted  at  high  blow  parameters  and  we  failed  to  obtain  a  steady 
solution  for  this  level  of  regime  parameters  (jw  I  PqUq  ~  0,0l) .  A  similar  problem  was 
discovered  in  paper  [14],  and  calculation  was  possible  only  with  use  of  the  model  of  mixing 
length. 

In  order  to  compare  the  results  with  experimental  data  solution  of  incompressible  gas 
dynamics  equations  were  performed  without  boundary  layer  approximation.  Boundary 
conditions  were  taken  from  experimental  data  described  in  [5].  Simulation  and  experimental 
results  are  shown  in  Fig.  6.  From  Fig.  6a  we  can  see  that  simulation  gives  correct  results  on  the 
flame  position  while  temperature  in  the  flame  front  is  20%  lower  than  in  the  experiment.  Fig. 
6b  and  Fig.  6c  show  good  agreement  between  simulation  and  experimental  data  on  mean 
longitudinal  velocity  and  its  pulsation. 
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